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Abstract. I review the basic observational properties of accreting millisecond pulsars that are 
important for understanding the physics involved in formation of their pulse profiles. I then discuss 
main effects responsible for shaping these profiles. Some analytical results that help to understand 
the results of simulations are presented. Constraints on the pulsar geometry and the neutron star 
equation of state obtained from the analysis of the pulse profiles are discussed. 
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INTRODUCTION 

First evidences for rapid rotation of neutron stars (NS) in low-mass X-ray binaries 
appeared in 1996 with the discoveries of sub-kHz coherent oscillations observed during 
X-ray bursts by the Rossi X-ray Timing Explorer (RXTE) [see Is, for review]. Two years 
later, the first accreting X-ray millisecond pulsar (AMP) was discovered [2], which was 
a nice confirmation of the theory of production of millisecond radio pulsar in the course 
of accretion. By 2008, there are at least 8 accreting NS discovered that show coherent 
oscillations for the extended periods of time during the outbursts [Hi Hi- These sources 
are exceptional laboratories to study the physics of accretion onto magnetized stars. 
They also have a great potential to test the NS structure. These goals can be achieved 
by studying the pulse profiles. This, however, requires understanding of the processes 
responsible for the production of the X-rays at the NS surface as well as detailed 
modeling of the propagation of the radiation to the observer. 

I review here the main observational data that have direct relation to this topic. These 
include averaged and phase-resolved broad-band spectra, pulse profiles and their energy 
dependence, and time lags. Then I describe recent efforts to develop theoretical models 
devoted to modeling these data and the main results obtained from these studies. 



OBSERVATIONS 

Broad-band spectra 

In order to understand how the X-ray pulses are produced in AMP, we need first 
to identify the basic radiation processes responsible for the X-rays we observe. The 
broad-band coverage of the RXTE together with XMM-Newton and INTEGRAL gives a 



possibility to study the spectra of AMP in great details. The AMP turn out to have rather 
similar, power-law like spectra with photon index F ~ 1.8 — 2.1 with a cutoff around 100 



keV 10, la B IsL B H \lM (see Fig. |t])lJ Such spectral shape strongly argues in favor of 
thermal Comptonization (with electron temperature kTg ~ 50 keV and Thomson optical 
depth of Tt ~ 1) as the main emission process, while the absence of strong spectral 
variability can be used as an argument that the emission region geometry does not 
depend much on the accretion rate and the magnetic field strength. The seed photons for 
Comptonization are probably produced by reprocessing of the hard radiation at the NS 
surface (similarly to the two-phase model developed for Seyfert galaxies [13, [3 (3), 
and not cyclo- synchrotron radiation. Below ~5 keV, a soft excess is often present. It 
can be associated with the blackbody like emission (with kT ~ 0.4 — 0.7 keV) from the 
heated NS surface, but not the accretion disk, because it is pulsating. 
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FIGURE 1. Broad-band spectra of SAX J1808.4-3658 (left) and XTE J1751-305 (right). Dotted and 
dashed curves show the Comptonized and blackbody components, respectively. The reflection is shown 
by triple-dot-dashed curve and an additional soft disk component is shown by dot-dashed curve. 

It is also of interest to study spectral variability as a function of pulsar phase. The 
simplest model has only two parameters corresponding to the normalizations of the two 
components (blackbody -i- Comptonization). One sees a hysteresis-like loop (see Fig.|2l), 
with the blackbody lagging the Comptonized emission laS]. This behavior results in 
the soft time-lags (discussed below). The oscillation amplitude and the pulse shape of 
the two components are also significantly different. 



Pulse profiles 

Most of the pulsars show extremely simple pulse profiles, which can well be described 
by just one sine wave. However, at high energies above 10 keV contribution of the 
harmonic becomes stronger (with the exception of IGR J0029 1-1-5934, [fiol]'). the profile 
becomes more skewed with faster rise and slower decay. During the decay of the 2002 
outburst of SAX J1808.4-3658, the profile show double peaks, which is probably the 



Spectral fits of the Chandra or XMM-Newton data alone in a narrow energy band give much harder 
slopes with T ~1.4-1.5 [see e.g.[T2l[T3ll. 
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FIGURE 2. Pulsation of the two model components, blackbody and Comptonization, in SAX J1808.4- 
3658; from 1^. Similar behavior is seen in XTE J1751-305, see |i8J. 
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FIGURE 3. Pulse profiles observed from SAX J 1808.4-3658 during its 2002 outburst and the best-fit 
sinewaves. (From Ibragimov et al., in preparation.) 



evidence for appearance of the secondary spot. Strong profile clianges also result in a 
very noisy behavior of the pulsar phase jlTl llSll. 



Time lags 



Pulse profiles at higher energies reach their peaks at an earlier phase relative to 
the soft photons resulting in the soft time lags jlol . 0, [sl, [Zll], with the absolute 
value increasing with energy up to about 7-10 keV and saturating afterwards. Only in 
IGR J0029 1+5934, the trend seem to reverse above 10 keV ['22!/l0]. Here it is important 
to note that the blackbody's contribution decreases exponentially with energy and around 
7 keV it is negligible. Thus, changes in the pulse profile of the two components are the 
obvious cause of the lags. One cannot separate lags from the pulse profiles and any 
successful time-lag model should actually explain not the lags, but the difference in the 
pulse profiles of the spectral components. A difference in the emission pattern [7] can 
be responsible for that. On the other hand, models involving Compton down- (or up-) 
scattering [|19l. l23|] pay attention only to the lags ignoring the pulse profiles H 



THEORY OF PULSE PROFILES 

The observed pulse profiles depend on a number of parameters: 

• Geometrical parameters, inclination / and co-latitude G of the magnetic poles. 

• Geometry of the accretion shock and the hot spot at the NS surface, which in turn 
depend on the accretion rate M, magnetic field strength and . 

• Emissivity pattern of the radiation, which may depend on M. 

• Stellar spin v, which affects the Doppler boosting and aberration effects. The 
importance of time delays, related to the finite propagation time around a NS, also 
depend on v. 

• Stellar compactness, i.e. NS radius-to-mass ratio R/M, which has a dramatic effect 
on the strength of gravitational bending. 

• Stellar oblateness, which might become important for stars of large radii rotating at 
V > 600 Hz. 

• Accretion disk inner radius, which blocks radiation from the antipodal spot and 
also reflects hard radiation. 

• Accretion stream's optical depth, which depends on M. The stream can eclipse the 
hotspot at some pulsar phases. 

The quality of the data does not justify yet to account for all possible effects. Many 
models make serious simplifications, for example, considering circular spots, spherical 
stars, blackbody emission, etc. 



^ A specific time-lag model of f23] suffers also from other problems. There the soft lags are explained by 
multiply scattering of the hotspot's hard radiation in the accretion disk. However, the radiation above 10 
keV is actually mostly reflected, while the few keV photons are immediately absorbed as the disk is rather 
cool. Thus, the role of multiply scattering is negligible and the lags cannot possibly be produced this way. 



Main effects 



Let us now briefly describe the main effects shaping the pulse profiles of AMP. 
Assume a small blackbody spot. Without general or special relativistic effects, a spot 
would produce sinusoidal variations (with possible eclipses) due to a change of the pro- 
jected visible area. The gravitational light bending tends to decrease the pulse amplitude 
[24, 25, 26, 27, 28, 29] (compare dotted and dashed curves in Fig. S]), while the pulse 
remains almost sinusoidal [129|] . The importance of Doppler boosting and light travel 
time delays was pointed out in IboIi . who also showed that the frame dragging affects 
the pulse profiles only at a ~ 1% level. Therefore, one can use Schwarzschild metric. 
To account for the Doppler effect one makes Lorentz transformation from the frame ro- 
tating with the star to a non-rotating frame and then follows the light trajectory in the 
Schwarzschild space-time to infinity ("Schwarzschild -i- Doppler" or SD approach). The 
full formalism for computing the pulse profiles accounting for Doppler boosting, rela- 
tivistic aberration and gravitational bending appeared for the first time in |0] (see [31] 
for detailed derivations). Accuracy of the SD approach is discussed in jjsll . Comparing 
to the slowly rotating stars, the projected area is changed by the Doppler factor 5 due to 
aberration, and the bolometric intensity is multiplied by 5^ [see 7, 31]. Since 5 reaches 
the maximum a quarter of the period earlier than the projected area, the pulse becomes 
skewed to the left (see Fig.S]). Light travel time delays slightly modify the profile further. 
Oscillation waveforms and amplitudes were computed in SD approach in [|32l.l33l]. who 
however present no formulae. The graphs shown in [[32L [33l] can be understood easily 
using analytical formulae for oscillation amplitudes derived in Hill (see Fig. [51). 

The light curves from realistic spots produced by accretion onto inclined magnetic 
dipole are computed in [|36l] using formalism o f A lso stellar oblateness might 
affect the pulse profiles for rapidly rotating stars 1371. l38]. Effects of the anisotropy of 
Comptonized radiation from the shock on the light curves and polarization were studied 



Angular pattern 

The angular pattern of the radiation escaping the NS surface plays an important role 
in shaping the pulse profile. The shape of the broad-band spectra of AMP implies that 
Comptonization happens in the optically thin atmosphere. So far only static atmospheres 
of constant temperatures have been considered [Sll^- As blackbody photons from the 
NS surface are scattered in the hot shock, they gain energy. The angular distribution 
of photons strongly depends on the number of scattering (and therefore their energy), 
saturating after about 5 scatterings. Radiative transfer calculations show that the angular 
dependence of the specific intensity can be approximated by a linear relation /(/i) oc 
1 + a/i (where ji = cos a and a is the angle between photon propagation direction 
and the local normal). For optically thin shocks, the scattered radiation correspond to 
negative a < —0.5. This gives pulse profiles dramatically different from those obtained 
assuming the blackbody (i.e. a = 0) or the Hopf profile (i.e. a — 2) corresponding to the 
optically thick electron scattering atmosphere [j4^ l4l|]. Changes in the emission pattern 
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FIGURE 4. Pulse profile for a small blackbody spot at the neutron star surface, (a) Dotted curve 
shows the profile for a slowly rotating star in Newtonian approximation. Gravitational light bending 
reduces variability amplitude (dashed curve). Doppler effect due to rapid rotation skews the profile 
(dot-dashed curves). Accounting for the light travel delays further modifies the profile slightly (solid 
curve), (b) Doppler factor as a function of phase. We assumed neutron star mass M — 1.4M0 and radius 
R = 2.5Rs = 10.3 km (Rs = 2GM/c^), rotational frequency V = 600 Hz, inclination ; = 45°, and the 
colatitude of the spot center = 45°. From 101 • 



happen smoothly with energy until blackbody flux becomes negligible (above 7-10 keV, 
see Fig. [l]) and only scattered radiation is left. Correspondently, the pulse profile shape 
and the position of their maxima moves to earlier phase [7]. This results in the soft time 
lags as observed. Thus Comptonization model explains both the pulse profiles and the 
time lags. I note that the lags due to finite travel time in the Comptonizing medium are 
completely irrelevant here. 



As the pulse profiles depend in a complicated way on the number of parameters, it 
is easier to understand their behavior using analytical approximations. In an important 
paper, Beloborodov [29] showed that in Schwarzschild metric there is a simple relation 
between the angle the photon momentum makes with the normal a and that angle 
without the light bending i// [see alsol28ll: 



where u = Rg/R, cos i// = cosz cos + sin z sin cos ^, and ^ is the pulsar phase. With 
the same accuracy the time delays can be computed using a relation A? = ( 1 — cos i//)7?/c 
and the Doppler factor can also be approximated as [31] 



Analytical approximations 



cos a ^ M + ( 1 — m) cos V/, 



(1) 



5 = 1 — /3eqVl — u sinz sin COS 0, 



(2) 
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FIGURE 5. Amplitudes of the fundamental and the harmonic for the pulsars with two antipodal spots 
as a function of magnetic colatitude 9. Note that for / = 90°, only a harmonic is present. Solid curves are 
exact calculations and dashed curves are analytical approximation for M ~ 1.4M0, R = 2.5Rs, V = 600 
Hz, and a = -0.7. Adapted from Il3ll1 . 



where /3eq = Veq/c = InRv/cy/l — u is the equatorial velocity. This allows to compute 
the pulse profiles, Fourier amplitudes and phases analytically with high accuracy (see 

Analytical formulae allow to understand easily that for one spot the pulse amplitude 
increases linearly with sin? sin0. The deviation of the pulse profile from the sine- wave 
(i.e. the ratio of the harmonic to fundamental) also depends linearly on sinz sin0, the 
rotational velocity /3eq ^ 1 as well as on the anisotropy parameter a (which does not 
need to be small). Thus strong harmonics (skewness of the profile) are mostly a result of 
the emission anisotropy. We expect strong harmonics only if the pulse itself has a large 
amplitude (as e.g. in XTE J1814-338, [42]). 

It is known that, for the blackbody emitting spots, oscillations almost disappear when 
both spots are visible (for the so called class IV pulsars, [|29|] '). because each spot 
produces a sine wave which are shifted in phase by 180°. An anisotropic source with 
a 7^ 0, however, produces strong oscillations. The amplitudes of the fundamental and 
harmonic behave roughly as oc a sin2/ sin20 and oc a sin^z sin^0, respectively i^M, 
with a weak dependence on stellar compactness and rotational frequency. It is clear thus 
that even a slight anisotropy is capable of producing strong pulsations. Stellar rotation 
alone, even for a = 0, also produces oscillations with amplitudes oc sin2z sin 20 and 

oc /3eq sin^ / sin^ 0. This argues in favor of weak magnetic field (and not their large mass 
as proposed recently in ll43ll ) in the majority of NS in low-mass X-ray binaries as the 
main reason for rarity of AMP. Alternatively, a nearly perfect alignment of the magnetic 
and rotational axes, ~ 0, solves the problem [|44|1 . 
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FIGURE 6. Constraints on the inclination and magnetic inclination for XTE J1751-305. The lines 
; = 30° and i = 85° are the lower and upper limits on the inclination of the system ll45ll . The absence 
of eclipses of the primary spot limits the allowed region to cos(/ +9) > —K (where K = R/{R — R^); to 
the left from the corresponding line, see [29J for description of the regions). The antipodal spot can be 
invisible if it is blocked by the star, this gives the allowed regions of parameters i and 9, cos(/ + 9) > K 
(lower hatched region). It can also be blocked by the accretion disk (inner disk radius of Ri„ = 40 km 
is assumed), this gives another allowed region (upper hatched area). Contours of constant variability 
amplitude A — 0.045 are shown for the blackbody spot at a slowly rotating star (dashed curve) and for the 
Comptonized radiation (dotted curve) multiply scattered in the shock of Thomson optical depth Tt = 1 .5 
and a star rotating at v = 435 Hz (as in XTE J1751-305). NS of mass M = 1.4M0 and radius R = 3Rs is 
assumed. From jsl]. 



Analytical approach gives the peak-to-peak amplitude (for a> —1/2, small spot and 
slow rotation) of nearly sinusoidal pulses lt3lll 

(3) 



^max ^min 

U{l+2aQ) 

+ 



where Q = u + {I — u) cos /cos 0,U = (1 — m) sinz sin0. This simple relation allows us to 
understand the general behavior of the oscillation amplitude even for the rapid rotation, 
which can be used to put constraints on i and 6 (see Fig. (6]). 



Accretion disk: eclipses and reflection 

As the inner edge of the accretion disk in AMP reaches a few stellar radii, the 
secondary accretion spot is invisible for almost any inclination. This is supported by 
the nearly sinusoidal pulse profiles observed in most AMP. At low accretion rate, the 
Alfven radius grows and the disk retreats towards the corotation radius allowing direct 
view of the antipodal spot. The profile becomes much more complicated (as seen in 
SAX J1808.4-3658, Fig. [S]), which results in the phase jumps of the fitted sinusoids. 
This possibly explains a complicated behavior of the pulsar phase [ iTlflsl] . 



Using simple analytical formulae for light bending derived in (t^, one can get 
constraints on the inclination and magnetic angle from the absence of the secondary 
spot [8]. It can be blocked either by the accretion disk implying rather high inclination 
(see upper hatched region in Fig. (6]) or by the star itself (lower hatched region). However, 
if double peaked profiles are sometimes seen (as is the case for SAX Jl 808.4-3658 at 
low luminosity, see Fig. [3]), only the first solution is possible. 

The accretion disk not only absorbs some radiation, but also reflects hard X-rays. The 
reflection amplitude depends on the inner disk radius and stellar compactness ^ (as 
well as on the emission pattern and spot size and position). Variation of the reflection 
with the accretion rate can be used to track the position of the inner disk radius. 

Constraints on the NS equation of state 

Detailed fitting of the AMP pulse profiles allows us to constrain the NS equation 
of state. Analysis of the data on SAX J 1808. 4-365 8 resulted in constraints on the NS 
radius S < R < 11 km for M = 1.4-1.6M0 [j7|]. Recently, these results were confirmed 
with the model including also time delays and stellar oblateness ^46"]. However, a study 
of XTE J1814-338 [47] gave a large NS radius requiring a stiff equation of state. This 
contradiction possibly originates either in the assumption of the spot geometry and/or 
simplicity of the emission model. Further studies on how these assumptions affect the 
results are required. 

CONCLUSIONS 

I have reviewed the basic observational properties of AMP such as spectral energy 
distribution, pulse profiles, time lags that are important for understanding of the physics 
of accretion in these sources. I also discussed the major effects that are responsible 
for shaping the pulses in AMP: light bending, Doppler boosting, anisotropy of emission, 
eclipses by the accretion disk, etc. The complicated dependences of the pulses on various 
parameters trigger the analytical approach which allows to understand better the results 
of simulations. I discussed how various observables such as amplitudes of fundamental 
and harmonic depend on the main parameters as well as how the invisibility of the 
secondary spot constrains geometrical parameters and the inner disk radius. This in 
turn can be used to constrain the pulsar magnetic field. Various attempts to measure 
NS compactness produce contradicting results, which probably reflects our ignorance of 
the angular distribution of the radiation as well as the spot geometry. Theoretical studies 
of these aspects of accretion onto AMP are certainly welcome. 
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